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Influence of AGN jets on the magnetized ICM 
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ABSTRACT 

Galaxy clusters are the largest structures for which there is observational evidence of 
a magnetised medium. Central cores seem to host strong magnetic fields ranging from 
a few 0.1 /J,G up to several 10 /xG in cooling flow clusters. Numerous clusters harbor 
central powerful AGN which are thought to prevent cooling flows in some clusters. 
The influence of such feedback on the magnetic field remains unclear: docs the AGN- 
induced turbulence compensate the loss of magnetic amplification within a cool core? 
And how is this turbulence sustained over several Gyr? Using high resolution magneto- 
hydrodynamical simulations of the self-regulation of a radiative cooling cluster, we 
study for the first time the evolution of the magnetic field within the central core 
in the presence of a powerful AGN jet. It appears that the jct-induccd turbulence 
strongly amplifies the magnetic amplitude in the core beyond the degree to which it 
would be amplified by pure compression in the gravitational field of the cluster. The 
AGN produces a non-cooling core and increases the magnetic field amplitude in good 
agreement with /.jG field observations. 

Key vifords: galaxies: clusters: general ~ galaxies: cooling flows - galaxies: jets - 
galaxies: magnetic flelds - methods: numerical 



1 INTRODUCTION 

There is increasing evidence for the strong magnetization of 
the hot plasma of the i ntra-cluster medium (ICM, see the 
review of lGovoni fc Ferctti 2004 1 , the largest scale magnetic 
field that has been constrained. But the evolution and the 
origin of cosmological magnetic fields remain open questions. 
Cosmological magnetic fields could have a primordial origin 
or could be propagated by galactic winds emerging from 
galaxies with strong dynamos. The lack of observational in- 
formation on the magnetisation of the inter-galactic medium 
keeps us from answering these questions. 

Hopefully, numerical simulations may shed light on 
these issues. It is starting to become clear that mag- 
netic fields in galaxy clusters are compressed essentially 
by the gravitational collapse of the gas, and are sub- 
stantially amplified by the shear motions developed in 
this hot phase by the turbulent motions of the shock- 
heat ed surfaces, the galaxy motions and minor / major merg- 
ers (iRoettiger et al.j [l999l: lOolag et al.l Il999l, |2002|. |2005|; 
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Clarke et al.l 2001 ) and cool cores (iTavlor fc PerlevI 



Vogt fc EnfilinI l2003l : lEnfilin fc Vogtl |2006| ): cooling 
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.Subramanian et alT 20061 : Asai et ahl 

20071 : iDubois fc Tevssieij|2008t ). Mean magnetic fields inside 
cluster cores range from a few 0.1 /iG up to several 10 fiG 
for the most massive ones, refiecting the natural scaling 
between magnetic field and density. There is al so a difler- 
ence between magnetic fields in non-cool cores (|Kim et al.l 
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flows seem to be a d river of the magnetic amplification 
ICariUi fc Tavloij|2002l ). 

Cosmological numerical simulations are a ble to repro- 
duce these trends, adiabatic simulations from iDolag et al.l 
(|2005l ) showing that the largest clusters are permeated with 
the strongest mag n etic fi elds. Radiative simulations from 
[Dubois fc Teyssien (|2008l ') have demonstrated that a cool- 
ing fiow within a galaxy cluster leads to a dense cool core 
and increases the quantity of turbulence and magnetisation 
in the inner 100 kpc, thereby explaining the fundamental 
difference between cool cores and non-cool cores. But the 
absence of cooling flows in some clusters is still puzzling as 
all cores should have already experienced catastrophic col- 
lapse jPabian. 199-1 ). 

There are several possible explanations of this puzzle: 
central cores could have b een heated by therma l conduc- 
tion from the outer parts (|Voigt fc FabianI 12004), or they 
could have been pre-heated by stellar outbursts ([Babul et al.l 
|2002| ). But the most popular explanation is that a cen- 
tral pow erful Active Galactic Nuclei (AGN) is re-hea ting 
the core (|Binnev fc Taboilll995l : [kephaeh fc Silklll995l ). X- 
ray surveys show remarka bly high energetic outfiows in the 
centr e s of galaxy cluster s jArnaud et al. 19841: Carilli et all 
1994 iMcNamara et ai] 1200 ll. |2005|; iFabian et al.l l2o'oa 



Birzan et a l. 2004; Forman et al.ll2007^ . Such bursts of en- 



ergy are associated with supermassive black holes (SMBHs) 
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accreting hot gas and propelling sup ersonic jets into the am- 
bient medium of the ICM plasma (|Progall2003l : iMcKinnevI 
I2OO6I ). 

Different attempts have been made to determine the ef- 
fect of central feedback on the ICM jBriiggen fc Kaiseij|2002l: 



Heinz et al]|2006l'l and on SMBH growth (jPi Matteo et aP 



2005l:ISiiacki et al in clusters hosting a powerful AGN 



Cattaneo fc Tevssieij (|2007h have pointed out that accretion 



onto the central SMBH is a self-regulated process and that 
the AGN prevents the formation of a cool core: the jet en- 
ergy release in the ICM is sufficient enough to re-heat the 
ICM. 

In this Letter, we address the question of the evolution 
of the magnetic field in the presence of a self- regulated AGN. 
The exact role of this feedback on cluster core magnetic 
fields is still ambiguous. Heating the core might decrease 
the central density and thus the magnetic energy, exerting 
a negative feedback on the field, but on the other hand the 
jet could increase the turbulence and non-negligibly amplify 
the overall central magnetic field, thereby exerting positive 
feedback. 

This problem is studied for the first time in an idealised 
(but cosmological) context, providing new insights on mag- 
netic field evolution inside galaxy clusters. 



2 NUMERICAL CONSIDERATIONS 

We are using t he Adapt i ve M esh Refinement (AMR) 
code RAMSES (|TevssieJ I2OO2I ') to solve the full set 
of ideal magneto-hydrodynamics ( MHD) equations us- 
ing a constrained t ransport method (|Tevssier et al.l I2OO6I : 
iFromang et~al]|2006l ') to prese rve the divergence of the mag- 
netic field. The HLLD solver (jMivoshi fc Kusandlioosi ') used 
in this simulation gives a five wave MHD solution, neglecting 
the contribution of the two slow magneto-acoustic modes. 

The initial conditi o ns ar e the same as the ones used 
in lCattaneo fc Tevssieij (|2007l ). therefore for full details the 
reader is invited to refer directly to the latter paper. Let 
us revi ew the basics. Dark matter is modelled with a static 
NFW (|Navarro et al.|[l996l ) profile with a concentration pa- 
rameter c = 5.53 and a virial mass Mvir = 1.5 X 10" Mq. 
The gas density profile is computed assuming that the gas is 
in hydrostatic equilibrium with a baryon fraction /b = 0.1, 
and that it follows a polytropic equation of state with a 
polytropic index 7p = 1.14. We negl ect the cluster's rot a- 
tion, which in any case would be low (|Bullock et al.ll200lf ). 

We assume that the initial configuration of the magnetic 
field is the result of the pure collapse of a gaseous sphere, 
and that any other amplification process is negligible. This 
means that the magnetic amplitude must follow B oc p^^'^. 
To ensure that the divergence of the field is preserved V.B = 
0, we require the magnetic vector A to satisfy V.(V x A) = 
where 



A{xi,yi,Zi) oc 



p2/3 



{Xi, Zi) 





(1) 



at the centre {xi,yi,Zi) of the edge of each cell over the 
computational domain. This particular choice of scaling the 
magnetic configuration to p^^^ is justified by the fact that 
the magnetic field within the central 200 kpc closely follows 



the spherical adiabatic collapse (|Dubois fc Tevssied I2OO8I ). 
For this study, it is unimportant if the magnetic field in the 
outer regions is amplified by shear motions, because the jet 
has a limited range (as we will see in the next section) and 
only the central 200 kpc regio n is able to radiatively collapse 
(see [Dubois fc Tevssieij I2OO8I ') . We must point out that our 
results are independent of the adopted profile (here p^^^) 
as long as field saturation is not reached (~ 50 pG for a 
2.10^ K cluster). According to this previous paper, we are 
fitting the magnetic amplitude in the centre of the cluster 
to 10~^ fiG. The cosmological magnetic amplitude in the in- 
tergalactic medium is set to ~ 1.3 x 10~^ ^G at p —< p > 
i.e. r ~ 4 Mpc), for a fair compa rison with iDolag et al.l 
20051 ) and lDubois fc Tevssieil (120081 ) that are both using a 
10"^ ^G field in the ICM. 

Following the jet 



implementation 



ICattaneo fc Tevssieij (|2007j ). the SMBH is atest par- 
ticle of MsMBH = 3 X lO'' M0 accreting at a iBondil (|l952l ) 
rate Msmbh = 47r(GMsMBH)^p/Cs (whic h is a good approx - 
imation for X-ray emitting sources, see I Allen et al.ll2006l ). 
assuming that the sound speed Ca ^ where Vff is the 
infall velocity onto the SMBH. In this model, the SMBH 
releases its energy in purely mechanical jet form whenever 
AfsMBH ^ MEdd, where MEdd is the Eddington accretion 
rate. In this simulation, it is reasonable to follow this 
regi me since Msmbh < 0.2 Mr /yr and MEdd ~ 65 Mq /yr 
fsee ICattaneo fc Tevssieil l2007h . As we do not resolve the 
propagation of the jet, we assume that a certain quantity of 
the surrounding gas is entrained with the jet material as it 
propagates. We therefore increase the mass injection rate of 
the jet Mj = x Msmbh, by a mass loading factor given 
by 7; = 100. The jet releases e — 10% of the overall accretion 
energy in kinetic form such that the momentum rate of the 
jet is qj = -\/2eMsMBHC and its rate of energy deposition 
is E] = eMsmbhc^. The velocity of the jet can therefore 
be expressed as uj — qj/Mj — c\/2e/r] ~ 1350 km. s~^ 
in opposite directions. Mass, momentum and energy are 
spread over a small cylinder of radius rj = 3.2 kpc (5 cells) 
and height of 2/i,i where hj = 2.5 kpc (4 cells) multiplied 
with a kernel window function 



exp 



2r] 



(2) 



as m lOmma et all |200i). Height and radius are arbitrary 
chosen, but are sufficiently large such that the jet is sampled 
with several cells. Thus the jet is constrained to flow along 
the y axis with only an input Uy velocity component. 

The simulation takes into account the self-gravity 
of the gas within the dark matter potential. Gas is al- 
lowed to cool by radiating its energy via atomic collisions 
(jSutherland fc Dopitaj[l993l ) assuming 75 % Hydrogen and 
25 % Helium fractions. The hot plasma gas has a. Z = Zq/^ 
metallicity constant through time and space, allowing the 
gas to cool more efficiently down to 10 K. 

Our computational domain is Lbox =648 kpc long with 
a 64'' coarse grid resolution (^min ~ 6). We refine the grid 
according to a geometric strategy (no time evolution): the 
1 = 1 grid is embedded within a L]aax/9> sphere radius cen- 
tred on the SMBH, ^ = 8 and 1 = 9 within a Lbox/IB sphere 
radius, and I = 10 within a Lbox/32 sphere radius. In this 
way the resolution within the centre reaches an equivalent 
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Figure 1. Average magnetic amplitude (solid line) in the central 
core r < 50 kpc as a function of time for the no- AGN run. A pure 
p^^^ scaling law is ovorplottcd (dashed line). 

1024^ grid (^^ax = 10) allowing for a minimum cell size 
Ax = 0.64 kpc. 

3 RESULTS 

3.1 Cooling flow clusters 

We tackle first the case of no AGN in the galaxy cluster. 
Since there is nothing to prevent the core from collapsing, 
a strong cooling flow develops and the gas density reaches 
very high values. We compute the density-weighted average 
magnetic field as 

for each cell lying within the core radius taken to be rcoro = 
50 kpc. The definition of the core radius is quite ambiguous. 
However what is important here is to properly measure the 
excess of magnetic amplification relative to the case of pure 
collapse. As we weight the average magnetic amplitude by 
the density, denser regions contribute more to the results. 
Thus taking a larger core radius will not lead to a different 
behaviour, but care should be taken that the core radius is 
not so small that all the effects of the turbulent amplification 
(in the jet case) are lost. In the AGN case, we also checked 
that the magnetic field in the core is independent of the 
unresolved region of the jet by excluding this small region (< 

4 kpc) from the calculation of the density-weighted average 
magnetic field. 

As shown in fig. [T] the magnetic field rises up to a few 
^G in the no- AGN case, because of the catastrophic collapse 
occurring at 4 Gyr. As there is no turbulence in this run, the 
magnetic field rather nicely follows the evolution expected 
for pure compression (B oc p^^^ where p is the average den- 
sity). The discrepancy between the measured magnetic am- 
plitude and that predicted by pure compression is due to 
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Figure 2. Same as fig. [T] for the AGN run. 

spurious (numerical) magnetic reconnection in the centre of 
the core where the fiow converges. 

3.2 AGN feedback: non-cooling flow clusters 

If we make the same measurement for the AGN run, the 
results are fundamentally different (fig. the magnetic 
amplitude in the core reaches only 0.1 pG, but it is strongly 
amplified relative to the pure compression regime. Initially 
the magnetic field is a dipole-like structure aligned with the 
z axis, and the jet (aligned with the y axis) is perpendicular 
to the field. The jet stretches and amplifies the y component 
of the field in a few 100 Myr up to 0.1 pG, as one can see 
in fig. |3ja). The jet propagates far from the centre of the 
cluster (~ 100 kpc) until the ram-pressure of the ICM has 
sufficiently dissipated its mechanical energy. The cluster core 
and its magnetic evolution remains in this quasi-stationnary 
state. As there is no initial turbulence in the cluster, the jet 
remains nearly unperturbed. 

At ~ 6.5 Gyr, the cooling catastrophe occurs; the very 
centre of the cluster core collapses in a free-fall time. In this 
case, the jet is unable to prevent this critical collapse, but 
as the density rises in the centre, the accretion rate onto 
the SMBH grows and as a consequence the energy released 
by the jet increases up to the point that the jet begins to 
halt the accretion onto the SMBH. In this way, the clus- 
ter core evolution is self-regulated by its AGN. This self- 
regulated situation is highly unstable, and alters the prop- 
agation direction of the jet, leading to turbulent motions 
within the cluster core. Shear motions twist magnetic field 
lines (fig. 131(b)), so that in the core the magnetic field is 
no longer aligned with the initial field. The relatively strong 
magnetic amplitude on the z axis, at t = 8 Gyr and z > 100 
kpc, as opposed to that along the y and x axes is the result 
of the preferential orientation of the initial field (the same 
situation occurs if there is no AGN feedback) . 

We remark that the self-regulation process of the AGN 
is crucial for the magnetic field evolution. Each time the 
cluster core suffers a rapid collapse due to cooling (at t ~ 6.5 
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Figure 3. Projection along tlie x axis of the density-weighted average of the magnetic amplitude in fiG units after (a) t = 2 Gyr (left 
panel) and (b) t = 8 Gyr (right panel) for the AGN run. The image size is -I/box/4 = 162 kpc. In the left panel: magnetic field is stretched 
and amplified along the jet axis. In the right panel: catastrophic collapse has occurred, the jet propagation is no longer straightforward 
and brings turbulence into the plasma. Thus the magnetic field is filling a small sphere and some magnetic filaments are created due to 
shear motion of the gas. The vertical magnetic feature is due to the initial configuration of the field (along the z axis). 



Gyr and t ~ 8.2 Gyr), the overall magnetic amplitude in 
the core rapidly grows by compression. As the SMBH is fed 
at increasing rates, the AGN is strenghtened, leading to a 
violent outburst in the core which inhibits the cooling flow, 
powers the turbulence, and hence the magnetic amplification 
within the cluster core. 

All this amplification is mainly concentrated within the 
centre of the cluster as shown by fig. ID We can see that the 
magnetic field at t = 8 Gyr is strongly amplified (by a factor 
~ 4-5) up to a distance of ~ 100 kpc from the centre. As 
a source of turbulence, the jet can propagate shear motions 
only up to that distance, meaning that all the jet-driven 
turbulence is concentrated within the core. The typical 3D 
velocity dispersion in the core at t = 8 Gyr is g ~ 100 — 
170km/s in fair agreement with lSubramanian et al. I (|2006l ). 
Finally the AGN is sufficient enough to prevent the creation 
of a dense and cool core. 



4 DISCUSSION 

These numerical results of the magnetic evolution of a galaxy 
cluster with or without an AGN jet are able to explain the 
discrepancy between strong magnetic fields in cool cores (no 
AGN heating) and low magnetic fields in non-cool cores 
(AGN heating). Without any AGN, a strong cooling flow 
appears and the magnetic amplitude in the cool core (< 50 
kpc) reaches ~ 4 /iG by gravitational compression alone. The 
presence of a powerful AGN jet, on the other hand, leads to 
softer and hotter cores within which the magnetic ampli- 
tude is ~ 0.1 /iG but turbulence, rather than gravitational 
compression, drives the ampliflcation of the field. 

This work is in g ood agreement with previo us numerical 
studies. In particular. [Dubois fc Tevssied (|2008l ) have shown 
that adiabatic cores (no AGN feedback and no cooling) are 
unable to sufficiently amplify the magnetic field within the 
central core (10~^ fJ-G), but cooling flows are able to match 




r (kpc) 

Figure 4. Average magnetic amplitude (solid lines) and p^^^ law 
(dotted lines) at t = Gyr (thin lines) and t = 8 Gyr (thick lines) 
for the AGN run. One can see that the initial conditions for the 
magnetic field arc fitted on the p^^^ profile. 



observations (several 0.1 /iG). Here we have shown that AGN 
are able to stop a central cooling flow and to prevent the for- 
mation of a dense core, and at the same time the jet-driven 
turbulence can amplify inner magnetic fields up to a few 
nG in accordance with observations of non-cool cores. The 
scenario of a thermally conducting plasma in the absence 
of AGN feedback seems to be weakened. Thermal conduc- 
tion can re- heat the core, but it d oes not bring in sufficient 
turbulence (|Kim fc Naravanll2003l ). As we need strong tur- 
bulent amplification of the initial magnetic field for non-cool 
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core clusters, it is difficult to explain how fiG fields can be 
reached in scenarios that only rely on thermal conduction. 

Of course we are far from reaching the full magneto- 
turbulent cascade with the kiloparsec resolution of cosmo- 
logical simulations. But such simulations give access to the 
full evolution of a galaxy cluster over a fraction of the Hub- 
ble time, and allow one to explain the large-scale evolution 
of these magnetised structures. We must also point out that 
the only source of turbulence in our simulation is the jet 
itself, meaning we neglect all the contributions coming from 
the outer parts of the cluster. Shock-heated turbulence could 
be important, especially in destabilising the jet and forcing 
it to stay in the core. We also did not consider any magne- 
tization of the jet that could either coUimate it or magne- 
tize the entire cluster core. Indeed, even though the amount 
of magnetic e nergy inje c ted in to the ICM by AGN jets is 
yet unknown, IXu et ahl (|2009l ) have shown that a tempo- 
rary deposit of magnetic energy within the jet, coupled to 
an efficient turbulent dynamo is able to somewhat enrich the 
hot plasma. We simply claim in this letter that within the 
interior of clusters, self-regulated AGN jets are able to sus- 
tain turbulence in the long term and to therefore efficiently 
amplify/transport magnetic fields. 

This Letter is the first attempt to self-consistently con- 
sider the evolution of AGN feedback powering turbulence in 
a magnetised cluster. Our work suggests that it is of par- 
ticular importance to include AGN feedback if one wants to 
reproduce observations and explain the origin and the evolu- 
tion of the magnetic field on very large scales of the Universe. 
In a forthcoming paper, we will address this problem in a 
fully non-idealised cosmological context. 
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